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The Viking Gas Chromatograph Mass Spectrometer failed to 
detect organic compounds on Mars, and both the Viking Labeled 
Release and the Viking Gas Exchange experiments indicated a 
reactive soil surface. These results have led to the widespread belief 
that there are oxidants in the martian soil. Since H 2 0 2 is produced 
by photochemical processes in the atmosphere of Mars, and has 
been shown in the laboratory to reproduce closely the Viking LR 
results, it is a likely candidate for a martian soil oxidant. Here, 
we report on the results of a coupled soil /atmosphere transport 
model for H 2 0 2 on Mars. Upon diffusing into the soil, its concentra- 
tion is determined by the extent to which it is adsorbed and by 
the rate at which it is catalytically destroyed. An analytical model 
for calculating the distribution of H 2 0 2 in the martian atmosphere 
and soil is developed. The concentration of H 2 O z in the soil is 
shown to go to zero at a finite depth, a consequence of the nonlinear 
soil diffusion equation. The model is parameterized in terms of an 
unknown quantity, the lifetime of H 2 0 2 against heterogeneous 
catalytic destruction in the soil. Calculated concentrations are 
compared with a H 2 0 2 concentration of 30 nmoles/cm 1 * 3 , inferred 
from the Viking Labeled Release experiment. A significant result 
of this model is that for a wide range of H 2 0 2 lifetimes (up to 10 5 
years), the extinction depth was found to be less than 3 m. The 
maximum possible concentration in the top 4 cm is calculated to 
be —240 nmoles/cm 3 , achieved with lifetimes of greater than 1000 
years. Concentrations higher than 30 nmoles/cm 3 require lifetimes 
of greater than 4.3 terrestrial years. For a wide range of H 2 0 2 
lifetimes, it was found that the atmospheric concentration is only 
weakly coupled with soil loss processes. Losses to the soil become 
significant only when lifetimes are less than a few hours. If there 
are depths below which H 2 0 2 is not transported, it is plausible 
that organic compounds, protected from an oxidizing environment, 
may still exist. They would have been deposited by meteors, or be 
the organic remains of past life. <£> 1994 Academic Press, Inc. 

1. INTRODUCTION 

One of the most striking results from the Viking Lander 

missions to Mars was the lack of organic compounds 


found at the surface. The Gas Chromatograph Mass Spec- 
trometer (GCMS) showed that no organics were present 
in the martian soil at levels of parts per billion for complex 
organics, and parts per million for simple ones (Biemann 
el aL 1977, 1979). Even with the absence of life or the 
organic remains of life, this result is surprising. Organics 
should be supplied to the martian surface by meteoric 
infall, or could have been produced by an early reducing 
atmosphere and stored in subsurface reservoirs (Sagan et 
al. 1989). Recently, Flynn (1993) has estimated the flux 
of organic carbon on Mars to be approximately 5 x 10 18 
g cm -2 sec -1 . This is based on the assumption that 10% 
of the interplanetary dust flux is organic carbon and that 
carbonaceous material not subjected to a pyrolysis tem- 
perature of 900 K (Chyba et al. 1990) reaches the surface 
intact. The absence of detected organics near the surface 
implies the active destruction of these molecules, either 
photochemically or through oxidizing agents in the soil 
(Klein 1979). The oxidizing nature of the martian soil has 
also figured prominently in discussions of future explora- 
tions of Mars. The putative oxidants may be responsible 
for the destruction of organic material to considerable 
depth, thus precluding the recovery of reducing material 
that may be relics of early biological forms. Additionally, 
there have been serious expressions of concern regarding 
the effect that soil oxidants may have on materials used 
in robotic missions and on health and safety during human 
missions (Coulter 1987). The concern has centered on 
possible irritation of the respiratory system due to dust 
carried into the martian habitat through air locks. 

In this paper, we consider the production, transport, 
and destruction of an important possible martian soil oxi- 
dant, H 2 0 2 . A one-dimensional coupled soil-atmosphere 
model is developed to investigate the possible extent and 
concentration of H 2 0 2 in the martian environment. Data 
from the Viking biology experiments and from laboratory 
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TABLE I 

A Comparison of GEx 0 2 and LR 14 C Results 
(after Klein 1978) 


Sample 

GEx 0 2 
(nmol cm" 3 ) 

LR CO : 
(nmol cm -3 ) 

Viking 1 (surface) 

770 

—30 

Viking 2 (surface) 

194 

-30 

Viking 2 (subrock) 

70 

-30 


experiments that followed are used to constrain the model 
through appropriate boundary conditions. Our goal is to 
estimate the depth of the H 2 0 2 oxidizing layer, as well ns 
to point out the sensitivity of this result to the various 
assumptions and parameters used in the model. Further- 
more, we would like to show to what degree the soil is a 
sink for atmospheric H 2 0 2 , in order to assess the signifi- 
cance of surface interactions to atmospheric photochemi- 
cal models. 

The results of the Viking biology experiments have led 
to the widespread belief that there are oxidants in the 
martian soil. The key results upon which this hypothesis 
is based are (see Klein 1978, Mazur et al. 1978, Huguenin 
1982 for detailed reviews): 

• The GCMS failed to detect organics in surface samples 
and from samples below the surface (the maximum depth 
sampled was about 10 cm) (Biemann et al. 1977, Biemann, 
1979). Since there are at least two mechanisms that could 
produce organics on Mars, meteoritic infall and UV 
production (Biemann et al. 1977), the absence of 
organics suggests that a mechanism for destroying them 
is present. 

• The soil released 0 2 upon humidification in the Gas 
Exchange experiment (Oyama and Berdahl 1977, 1979) in 
amounts ranging from 70 to 770 nmoles cm' 3 (see Table 
I). The 0 2 was released rapidly within the first few hours 
after humidification; it then gradually tapered off and 
ceased. Heating of the sample to 145°C for 3.5 hr reduced 
the amount of 0 2 released by about 50% but did not elimi- 
nate it, indicating that biological activity was not involved. 
No additional 0 2 evolved when the samples were wetted 
with a nutrient solution but there was a slow evolution 
of C0 2 , indicating the oxidation of organics in the medium 
by an oxidizing agent. 

• In the labeled release experiment, i4 C0 2 was rapidly 
released when martian surface samples were wetted with 
an aqueous nutrient medium containing l4 C-labeled or- 
ganic compounds. The kinetics of the initial l4 C0 2 release 
were similar to the initial 0 2 release in the GEx humid 
mode tests. After the initial release tapered off (70 hr). 


14 C0 2 continued to evolve slowly until the experiment 
ended. The total abundance of C0 2 produced during the 
rapid response was suggestive of a nearly total reaction 
with only one component (formate) of the nutrient me- 
dium. A second injection of nutrient medium caused some 
30% of the I4 C0 2 gas to readsorb, followed again by its 
gradual release. The rapid ,4 C0 2 release was completely 
removed by heating at 160°C for 3 hr, partially destroyed 
at 40-60°C, and unaffected by storage for short periods 
at 18°C but lost after 2 to 4 months' storage at I8°C. 

H 2 0 2 is a good candidate for the thermally labile oxidant 
that produced the rapid evolution of C0 2 in the LR experi- 
ment. Hunten (1979, 1987) has suggested that H 2 0 2 is 
produced in the atmosphere by photochemical reaction 
at a rate of about 2 x 10 9 molecules cm" 2 sec" 1 and that 
this could be the source of the LR oxidants. Oyama and 
Berdahl (1979) duplicated the response of the LR experi- 
ment by reacting formate with H 2 0 2 and yFe 2 0 3 mixed 
with a suite of possible martian surface minerals. They 
also duplicated the slow C0 2 production in the LR and 
GEx wet cycle with a mixture of yFe 2 0, and formate. 
Ponnamperuma et al. (1977) also found that 14 C0 2 was 
produced when the Viking LR nutrient mixture was added 
to yFe 2 0 3 that had been irradiated with ultraviolet light 
under a Mars analog gas mixture. Heat pretreatment at 
160°C reduced but did not eliminate the effect. They sug- 
gested that H 2 0 2 may have been formed during irradiation 
which then reacted with formate in the medium. However, 
in other simulations. Levin and Straat (1981) found that 
H 2 0 2 reacted with other compounds in the nutrient me- 
dium besides formate. Thus the H 2 0 2 hypothesis does not 
account for the fact that only one component in the LR 
nutrient medium was oxidized to C0 2 in the initial fast 
reaction. Furthermore, Levin and Straat (1981) argue that 
H 2 0 2 is much more thermally labile than the oxidant that 
produced the LR result. 

In addition to photochemical formation of H 2 0 2 in the 
atmosphere, Huguenin et al. (1979) and Huguenin (1982) 
have suggested that chemisorbed H 2 0 2 is produced by 
frost weathering of olivine . Experiments with frost weath- 
ered olivine (Huguenin et al. 1979, Hugeunin 1982) repro- 
duced the release of 0 2 upon humidification (GEx result) 
and the oxidation of labeled formate solution to C0 2 (LR 
result), although the magnitude of the response was much 
larger than seen on Mars. Huguenin et al. (1979) argued 
that the similar kinetics of the release of C0 2 in the LR 
and 0 2 in the GEx experiments suggests that the same 
active agent is responsible for both reactions and pro- 
posed that the GEx reaction was caused by the catalytic 
dissociation of H 2 0 2 in the presence of iron catalyst and 
water. They further argued that the chemisorbed H 2 0 2 
would decay in the heat sterilization procedure to ad- 
sorbed 0 2 , which was desorbed upon humidification. 
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However, they did not demonstrate the heat sensitivity 
of these results experimentally. 

In order to estimate the extent of the oxidized layer, 
and to evaluate the possible hazard that soil oxidants 
might pose to future missions, a standard model of martian 
soil oxidants is needed. Such a model should have a clear 
mechanism for formation on Mars and should be consis- 
tent with the Viking biology experiment results. In this 
vein, H 2 0 2 has emerged as a “standard' 1 martian oxidant 
due to its formation in the atmosphere (Hunten 1979) and 
in the soil (Huguenin 1982) and its ability to reproduce the 
LR results, including their temperature sensitivity (Klein 
1978, 1979; Levin and Straat 1979). The concentration 
required to explain the LR results is about I ppm by mass. 
The arguments against H 2 0 2 as the soil oxidant include: 
(1) It has a short lifetime, 10 4 sec, against UV destruction 
on the martian surface. (2) It alone cannot explain the 
thermally stable GEx results. In the next section, we 
further develop the standard model and propose it as a 
useful benchmark pending more definitive information on 
the true nature of the soil oxidant. 

We investigate the diffusion and catalytic destruction 
of H 2 0 2 in the martian soil to determine whether its con- 
centration is consistent with the results of the Viking LR 
experiment. It is assumed that H 2 0 2 is produced in the 
atmosphere by photochemical reactions (Hunten 1979) 
and diffuses into the soil. The H 2 0 2 must be sufficiently 
mobile that it diffuses to depth before being photochemi- 
cally destroyed by surface UV light. We expect that the 
primary sink for H 2 0 2 is decomposition to water, probably 
catalyzed by the iron-rich soil grains. The depth of the 
oxidizing soil and the concentration of oxidant at the sur- 
face will depend on the degree to which H 2 0 2 is adsorbed 
onto soil grains and on the catalytic destruction rate of 
the H 2 0 2 . 

In the absence of the consideration of diffusion, adsorp- 
tion, and catalytic destruction of H 2 0 2 in the martian soil, 
atmospheric photochemical models (Hunten 1987, Kong 
and McElroy 1977) are only able to estimate a flux of 
H 2 0 2 into the soil. Kong and McElroy (1977) considered 
several ad hoc assumptions regarding the effectiveness 
of the soil as a sink for H 2 0 2 but did not attempt to quantify 
soil concentrations. Chyba et al. (1989) reported on a 
simple calculation of possible depths of H 2 0 2 in the mar- 
tian regolith, but did not consider the effect of adsorption 
and its feedback on deposition from the atmosphere and 
diffusion to depth. We concluded from these investiga- 
tions that in order to treat the problem of the depth to 
which H 2 0 2 may penetrate, a coupled soil-atmosphere 
transport model would be necessary. 

Presented here is an analytical model that predicts at- 
mospheric and soil concentrations of H 2 0 2 based on a net 
flux of H 2 0 2 into the soil and on adapting empirically 
derived adsorption isotherms for water on clean basalt 


under martian conditions (Fanale and Cannon 1971). The 
calculations are made for steady-state isothermal condi- 
tions at an average Mars surface temperature of 215 K. 
Concentrations in the atmosphere are determined from a 
simple mass balance equation and an exponential lapse 
rate. Concentrations in the soil are determined from a 
diffusion equation that incorporates catalytic destruction 
of adsorbed H 2 0 2 . Since diffusion is accomplished by 
the vapor phase, the adsorption isotherms are used to 
establish phase partitioning between the adsorbed (fixed) 
and vapor (mobile) phases. The soil diffusion equation is 
coupled to the atmoshpere equation through a flux bound- 
ary condition at the soil-atmosphere interface. The net 
flux of H 2 0 2 into the soil is established by the atmospheric 
concentration at the surface as well as by soil diffusion 
and loss processes. Similarly, atmospheric concentrations 
are dependent upon both atmospheric and soil processes. 

In order to uniquely determine a solution, the model 
requires that atmospheric production and loss rates, as 
well as soil loss and diffusion rates, be known. In addition, 
the equilibrium partitioning between H 2 0 2 phases mjst 
be known. The least well determined of these is the loss 
rate of adsorbed H 2 0 2 due to heterogenous catalysis. The 
lifetime of H 2 0 2 adsorbed in the soil is therefore taken as 
a free parameter, and the general solution with respect 
to this parameter may be calculated. 


2. THE MODEL 

The model is composed of two coupled partial differen- 
tial equations for the atmosphere and soil. Coupling is 
through a flux boundary condition at the soil-atmosphere 
interface, and must be consistent with atmospheric and 
soil concentrations. The phase partitioning between ad- 
sorbed and vapor phase H 2 0 2 enters in the soil equation. 


2.1 . Atmospheric Production 

Hydrogen peroxide is produced in the martian atmo- 
sphere through the combination of two hydroperoxyl rad- 
icals, 


H0 2 + ho 2 -+ h 2 o 2 + o 2 , 

with a rate constant k = 5.5 x 10" 12 cm 3 sec" 1 (DeMore 
1985). The production rate is calculated from 


</[H 2 Q 2 ] 

dt 


[H0 2 ][H0 2 ]*. 


Using a value of 6.5 x 10 7 cm 3 for the number density 
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of H0 2 and a scale height, H, of 18 km (Kong and McElroy The flux out of the soil is similarly calculated, but is 

proportional to n v (0), the vapor phase concentration of 
H 2 0 2 in the pore spaces near the surface: 


1977), the column production rate of H 2 0 2 is 
P = //[HO,] 2 /: = 4.18 x 10 K) cm" 2 sec' 


2.2. Atmospheric Loss 


F a = ^cn v (0). 


The loss of H 2 0 2 in the martian atmosphere is due Therefore, the net flux of H 2 0 2 into the soil is 
primarily to UV radiation at wavelengths below 3750 A, 
through photodissociation 


H 2 0 2 + hv 20H, 

with a rate coefficient (Shimazaki 1989) of 
/ = 2.5 x 10~ 5 sec 1 . 


F = F i ~ F a = t c« a (0) - ^t7j v (0). 


The net flux is represented mathematically through a dis- 
continuity in the number density, /i a (0) - « v (0), but it 
may also be related to the usual expression for diffusive 
transport, Fick’s Law, 


A simple calculation with these production and loss rates 
should yield an estimate of the atmospheric concentration 
at the surface, assuming a standard lapse rate, perfect 
mixing, and no losses to the soil. 


«.( 0 ) = jfi = 93 


10 8 


cm 


dN, 

dt 


= P -IN. d - F i + F 0 , 


(1) 


F = -D-. 
dz 


The abrupt change in number density may only be “felt” 
by molecules within a mean free path of the surface. 
Therefore, as far as the net transport of molecules across 
the surface is concerned, the effective concentration gra- 
dient is 


which is in rough agreement with the number density 
calculated by Kong and McElroy (1977), 2 x 10 9 cm 3 , 
using a more sophisticated photochemical model. 

2.3. The Atmosphere Equation 

If 7V a is the column density of H 2 0 2 in the martian 
atmosphere, 


dn 

dz 


« v ( 0) - nj 0) 


where X is the H 2 0 2 molecular mean free path. Fick’s 
Law therefore becomes 


D, 


F = -[n a (0) - n y ( 0)]. 


The diffusion coefficient, from kinetic theory, is 


where F l and F 0 are the fluxes into and out of the soil, 
respectively. The atmosphere is assumed to be isother- 
mal, with an exponential lapse rate and perfect mixing 


n A (z ) = n. d (0)e 


-ziH 


where « a (z) is the number density of H 2 0 2 as a function 
of altitude and H is the density scale height. Since we 
are interested in the long-term equilibrium distribution of 
H 2 0 2 Eq. (1) is solved in the steady state. The flux into 
the soil is assumed to be the rate at which molecules 
impinge upon the surface due to the thermal motion. If c 
is the average molecular speed, the flux into the soil is 


so that 


D -V k ' 


F = -c[nM - n v (0)]. 


where the difference in the multiplying factor is due to 
the fact that molecules approach the surface from all 
angles, and the appropriate spatial integration yields 1/4 
rather than 1/3. Solving for the atmospheric concentration 
at the surface, 


*a(0) = 


P + l/4cn v (0) 


F, = -cn a ( 0). 


!H + 1/4 c 


( 2 ) 
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2.4. Soil Diffusion 

It is assumed that the atmosphere is in pressure equilib- 
rium with gases in the interstitial spaces of the soil. With- 
out an external concentration gradient, the thermal flux 
of molecules impinging upon the plane of the martian 
surface would be equal to the thermal flux of molecules 
reentering the atmosphere. A concentration gradient is 
set up in the soil, however, by diffusion of the vapor 
phase downward, by the adsorption of H 2 0 2 onto grain 
surfaces, and by the subsequent catalytic destruction of 
the adsorbed phase. Diffusion of a gas through porous 
media may be modeled by two different mechanisms de- 
pending on whether the mean free path of the gas is greater 
or less than the average pore diameter. If pores are mostly 
larger than the mean free path, diffusion is dominated by 
molecular collisions. In this case, the diffusion coefficient 
is the gas-gas diffusion coefficient of H 2 0 2 in C0 2 
multiplied by a constant that depends upon the geometry 
of the soil. This transport mechanism is referred to as 
molecular diffusion. An analytical expression for the dif- 
fusion coefficient is obtained by modeling the pores in the 
soil as randomly twisted and distributed tubes embedded 
in a cross section of the soil (Ball 1981). Porosity is defined 
as the ratio of the pore space volume to total volume. 
The porosity of the martian soil is usually taken to be 0.5 
(Fanale el al. 1986). Tortuosity is the ratio of total tube 
length to length projected onto a preferred axis (the direc- 
tion of the flux). An average over an appropriate ensemble 
of tortuous tubes is used to characterize a given soil. The 
value used for most soils and thus for Mars is 5 (Satterfield 
1970, Smoluchowski 1967). The molecular diffusion coef- 
ficient is given analytically as D m = (e/t)D 0 where D 0 is 
the gas-gas diffusion coefficient, e is the soil porosity, 
and t is the soil tortuosity (Flood 1967). The value of 
the molecular diffusion coefficient for H 2 0 2 in C0 2 at an 
average surface temperature of 215 K is D m = 3.4 cm 2 
sec - 1 . 

If pores are mostly smaller than the mean free path, 
molecular collisions with the pore walls dominate the 
transport of H 2 0 2 , a process known as Knudsen diffusion. 
The average pore size may be very roughly taken to be 
about 1 pm (Fanale et al. 1986). The mean free path is 
calculated for a Mars average surface temperature of 215 
K. From the kinetic theory of ideal gases, it is found to 
be about 27.8 pm. Therefore, Knudsen diffusion is the 
dominant transport mechanism for H 2 0 2 through the mar- 
tian soil at this temperature. The Knudsen diffusion coef- 
ficient is given analytically as D K = (4(er 0 /T)[(2A:7Y 
7Ttn)] m , where r 0 is the average pore radius and m is the 
H 2 0 2 molecular mass (Clifford and Hillel 1983). The value 
of the Knudsen diffusion coefficient for H 2 0 2 in the mar- 
tian soil is calculated to be D K = 0.24 cm 2 sec -1 . The 
diffusion coefficient appropriate to the conditions consid- 
ered here may be derived by employing the Bosanquet 


approximation, where D 1 = D m l + D K l , and is found 
to be D = 0.22 cm 2 sec ] . 

2.5. Soil Loss 

One approach to determining the chemical lifetime of 
the soil oxidant is based upon the observed disappearar ce 
of the active agent in the LR experiment after 82 days of 
storage at 18°C. Levin and Straat ( 1979) used the observed 
decline in oxidant activity to infer activation energies of 
about 35-43 kcal mole 1 . This is much larger than typical 
activation energies for heterogeneous catalysis but is ap- 
proximately the activation energy (40-50 kcal mole ') or 
self-catalysis of H 2 0 2 at high temperatures (over 750 K; 
Schumb et al. (1955)). Chyba et al. (1989) have used tiis 
activation energy to compute a chemical lifetime for H 0 2 
in the martian soil of 10 s to 10 7 years. Laboratory work 
on the catalytic destruction rate of H 2 0 2 in a ferrous envi- 
ronment, however, indicates much shorter lifetimes, on 
the order of minutes (Schumb et al. 1955). Because of the 
uncertainty of the actual chemistry of the martian surface, 
catalytic destruction rates for H 2 0 2 under martian condi- 
tions are poorly constrained. The model is solved fo a 
wide range of lifetimes in order to assess the sensitivity 
of soil concentrations and atmospheric losses to differ- 
ences in soil H 2 () 2 lifetimes. 

2.6. Phase Partitioning 

The relative concentrations of vapor to adsorbed phase 
H 2 0 2 are calculated from an expression empirically de- 
rived by Fanale and Cannon (1971). The data was taken 
for the adsorption of H 2 0 onto clean basalt surfaces or 
concentrations several orders of magnitude higher than 
those expected for H 2 0 2 adsorbed onto martian soil 
grains. Nevertheless, since the dipole moment and heat 
of fusion per gram of H 2 0 2 are similar to those of H : 0 
(Weast 1974), we will use the Fanale and Cannon iso- 
therms as a starting point to estimate H 2 0 2 phase parti- 
tioning under martian conditions. The adsorption datu is 
represented by an expression that gives the number of 
grams of H 2 0 2 adsorbed per gram of soil as a function of 
temperature and partial pressure. The expression is 

p a = yP B *e~ m , 

where y, B 0 , and 6 are the empirically derived coefficie its 
with values y ~ 6.316 x 10 -9 , B 0 = 0.51, and 5 = 2679.8 
K. Using the ideal gas law and converting this to number 
densities, 

n c (z) = S(7>?», 

where n L (z) is the number density of the adsorbed phase. 
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n v (z) is the number density of the vapor phase, and 5(7) 
incorporates all the temperature dependence of the iso- 
therm data. In adapting these isotherms, we use the mass 
of H 2 0 2 rather than H 2 0, and use a value for the exponent 
of 0.5, rather than the measured value of 0.51. Then the 
phase partitioning function is 

nfz) = S(T)n ] J 2 , (3) 


2.7. The Soil Diffusion Equation 

For the soil diffusion equation, catalytic destruction of 
adsorbed H 2 0 2 is characterized by the lifetime r. The 
diffusion equation is then 

dn v (z, 0 dn c (z, t) 3F(z, t) n c (z, t) 

e : + — — = , (4) 

dt 6t dz T 


where, using p soi! = 1.3 g cm" 3 (Moore et al. 1977), at 
215 K, 


5(215) - 4.757 x 10 12 cm 3/2 . 

For a soil that is a poor atmospheric sink, the amount 
of adsorbed H 2 0 2 may be estimated, assuming that H 2 0 2 
in the soil pore spaces is in equilibrium with the atmo- 
sphere (/i v (0) — /i .,( 0 ) = 9.3 x 10 8 cm ' 3 ): \ 

n x = S(T)n l J 2 = 240 nmoles cm -3 . 


where e is the porosity of the soil and F is the flux. 
Transport in the soil proceed by diffusion, given by Fick’s 
Law, 


F(z, t) 


p dnfz, t) 
dz 


where D is the diffusion coefficient. Applying this to the 
diffusion term, and the phase partitioning to the loss term, 
the steady state soil diffusion equation is 


This represents a maximum amount of H 2 0 2 at the sur- 
face, assuming that soil losses are negligible. 

For the purpose of comparing the Viking LR experi- 
ment results with the model, the surface concentration is 
considered to be an average over the depth to which the 
martian soil was sampled. This depth is taken to be 4 cm 
(Klein et al. 1976). The maximum possible concentration 
of H 2 0 2 consistent with atmospheric-soil equilibrium was 
shown to be a factor of 10 higher than that measured by 
Viking. Therefore, soil diffusion and loss rates that do 
not deplete the surface by more than a factor of 10 will 
still be consistent with the Viking LR experiment results. 

In addition to solving the model in terms of the lifetime 
parameter, three cases are specifically considered. The 
model is solved in the Density Case by constraining the 
average surface concentration in the top 4 cm to be 30 
nmoles cm -3 . When the model is constrained in this way 
to be consistent with the Viking LR experiment results, 
we wish to see what lifetime is necessary to achieve such 
a steady-state surface concentration. We also wish to 
know how much of a sink the soil is for atmospheric 
H 2 0 2 under such constraints. The model is solved in the 
Lifetime Case for a lifetime of 10 5 years. By considering 
the adsorbed H 2 0 2 to be relatively stable on Mars, we 
would like to calculate the average surface concentration 
in the top 4 cm and again find to what degree the soil is 
a sink for H 2 0 2 . The third case examined is the Loss 
Regime Case. If H 2 0 2 is unstable enough, soil diffusion, 
driven by a strong gradient, and soil losses, could be 
sufficient to make the soil a significant sink for atmo- 
spheric H 2 0 2 . The lifetime required for the martian soil 
and atmosphere to be sinks of equal magnitude is calcu- 
lated. 


d 2 n, 5(7) l/2 
dz 2 Dr " v 


(5) 


with boundary conditions 

BC 1 lim F(z) = 0 and lim « v (z) = 0 

- large : large 


BC 2 


D 


dn v 

dz 


= Fi 


z = 0 


where the flux boundary condition couples to the atmo- 
sphere equation through 


Fj 



|l7! v (0). 


The nonlinear soil diffusion Eq. (5) may be integrated 
once, by multiplying by dnldz and assuming constant tem- 
perature, to give 


dn v 

dz 


45 
3 Dt‘ 


* 3/2 


+ Cl 


2.8. Boundary Condition I 

At some depth, we expect w v — > 0. We also expect 
the gradient to go to zero. Therefore, applying the first 
boundary condition implies that C, = 0. Then the concen- 
tration gradient becomes 
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dn v 
dz 




the general solution, Eq. (7). Other physically significant 
( 6 ) quantities may also be calculated from n c ( 0). They are 
presented in the following sections. 


It is expected that the concentration will decrease with 
depth. Therefore the gradient is negative. We may inte- 
grate again to obtain 


nfz) 



s y 2 i 4 

\2Dt) Z 


where C 2 may be identified as flv 4 ( 0 ) and must be deter- 
mined from the second boundary condition. Since we are 
interested in the concentration of adsorbed H 2 0 2 , we use 
the phase partitioning function, Eq. (3), to express the 
general solution to the soil diffusion equation in terms of 
the adsorbed phase concentration. 



2.10 . Extinction Depth 

The general solution to the soil diffusion Eq. (7), shows 
that the concentration of adsorbed H 2 0 2 goes to zero at 
a finite depth. Due to the nonlinear phase partitioning, 
n c = S(T)nl /2 , as the concentration of mobile H 2 0 2 de- 
creases with depth, the concentration of adsorbed H 2 0 2 
decreases at the square of this rate. The result is a rspid 
decrease in adsorbed H 2 0 2 with depth. If we denote the 
depth at which the concentration goes to zero as z e9 the 
extinction depth, then for n c (z e ) = 0 is given by, 


n' c /2 ( 0 ) - 


12Dr 


= 0 . 


So that the extinction depth is 


2.9. Boundary Condition 2 

From Eq. ( 6 ), the net flux of H 2 0 2 into the soil is, in 
terms of the second boundary condition. 


— D 


dn w 

dz 






in terms of atmospheric production and loss terms, 



< 4 (0) 


= P - 


IHn^i 0 ). 


Equation (2) for rt a (0), derived from the steady-state atmo- 
sphere equation is substituted into the above equation 



P - 


(l HP + (IH/4)cn v m 

IH + 7 <• 

4 



2.11. Average Surface Concentration 

The physically observable concentration of H 2 0 2 in the 
soil at the surface is the average of n c (z) over some depth. 
In the Viking Biology Experiments, the sample depth was 
approximately 4 cm (Klein et al . 1976). If the depth of 
the layer is Z\ and the average concentration in this layer 
is « h then 


1 r*i 

ni = - n Q (z)dz z x ^ £ e . 

Z\ Jo 

Since the concentration of adsorbed H 2 0 2 extinguishes at 
some depth, z e , the limits of this integration are valid cnly 
for Z\ < z e . For shallow extinction depths, where Z\ ^ z e , 
the average surface concentration becomes 


Again, since we are interested in concentration of ad- 
sorbed H 2 0 2 the phase partitioning function, Eq. (3), is 
used to convert to a polynomial in n c ( 0 ), 


« c 2 (0) + nfm 





PS 2 
IH 


= 0 , ( 8 ) 


which is solved numerically. 

Applying the second boundary condition has allowed 
us to calculate the adsorbed phase concentration for given 
lifetimes. The depth profile of H 2 0 2 in the soil is given by 


1 f *" e 

H| = - n Q (z)dz Z] 5: z e . 

Z\ J() 

Performing the integration, the average concentration of 
adsorbed H 2 Q 2 in the top 4 cm of soil is 
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2.12. Atmospheric Number Density 

Equation (2) for the atmospheric number density at the 
surface may be expressed in terms of the H 2 0 2 adsorbed 
phase concentration in the soil via the phase partitioning 
function, Eq. (3), 


«a(0) 


P + ~ 4 j,nl(0) 

IH + 7 c 
4 


( 11 ) 


Here, the atmospheric number density may be seen to 
explicitly depend upon both atmospheric and soil pro- 
cesses. 
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2.13. Net Flux into the Soil 

The net flux into the soil is given by the second bound- 
ary condition to the soil diffusion equation, 

= p, - p a = ( 3 ^)' ’ «v /4 <0). 

Once again, this quantity may be expressed in terms of 
the adsorbed phase concentration of H 2 0 2 through the 
phase partitioning function, Eq. (3), 



FIG. 1. The calculated depth profile of H 3 0 2 , with the average sur- 
face concentration in the top 4 cm constrained to be 30 nmoles/cm\ 
Due to the combined effects of diffusion and catalytic destruction of 
adsorbed H 2 0 2 , the concentration is seen to go to zero at a finite 
depth — approximately 1.5 cm. The model predicts an average lifetime 
of H 2 0 2 in the soil of 4.3 terrestrial years for this case. 


2 e — 1 .50 cm. 

Concentration as a function of depth for this case is plotted 
in Fig. 1. The atmospheric number density at the surface 
is calculated from Eq. (11) to be 

w a ( 0) = 9.170 x I0 8 cm" 3 . 


This completes the general solution to the model and 
allows us to compute the concentration depth profile, 
the average concentration in the soil at the surface, the 
atmospheric number density at the surface, and the net 
flux of H 2 0 2 into the soil, in terms of the soil lifetime 
parameter. 

3. RESULTS 

3.1 . Density Case 

If it is assumed that H 2 0 2 was responsible for the evolu- 
tion of C0 2 in the Viking Labeled Release experiment, a 
lower bound for the average concentration of H 2 0 2 in the 
top 4 cm may be taken to be about 30 nmoles cm -3 . 
This concentration is achieved in the model with a H 2 0 2 
lifetime in the soil of 

t = 1.35 x 10 8 sec = 4.3 terrestrial years, 

calculated from Eq. (10). The extinction depth is calcu- 
lated from Eq. (9) to be 


This may be compared with the number density that is 
obtained with the assumption that the soil is not a sink 
for H 2 0 2 , 

« a = 9.295 x 10 8 cm -3 . 

It can be seen that with an average surface concentration 
of 30 nmoles cm"\ the lifetime of 4.3 years yields a very 
shallow depth and a soil that is a poor sink for H 2 0 2 . 

3.2. Lifetime Case 

Assuming H 2 0 2 to be responsible for the Viking La- 
beled Release experiment results, Chyba et al. (1989) cal- 
culated the lifetime against catalytic destruction in the 
soil to be about 10^ years. Taking this as the constraint 
in the model, rather than average surface concentration, 
the same physical quantities are directly calculated. The 
resulting average concentration of H 2 0 2 in the top 4 cm 
is found to be 

n x = 236.6 nmoles cm -3 . 
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FIG. 2. The calculated depth profile of H : 0 2 , with the lifetime in 
the soil constrained to be 10 ? years. Although H 2 0 2 penetrates to greater 
depths than when the lifetime is 4.3 years, it can be seen that the 
extinction depth is still shallow — approximately 2.3 m. The average 
concentration of H 2 0 2 in the top 4 cm is calculated in this case to be 
236.6 nmoles/cm\ 


The maximum possible surface concentration was shown 
to be 240 nmoles cm -3 , so it can be seen that a lifetime 
of 10 5 years results in a surface concentration that is close 
to the asymptotic limit. The extinction depth is found to 
be somewhat deeper than for the Density Case, 

z e = 231.0 cm. 

H 2 0 2 concentration as a function of depth for the Lifetime 
Case is plotted in Fig. 2. The atmospheric number density 
reflects the fact that the soil is an even weaker sink for 

h 2 o 2 , 

/i a (0) = 9.29 x 10 8 nmoles cm“\ 

which is within 0.1% of that obtained for a completely 
passive surface. 

3.3 . Loss Regime Case 

The model results for the Density and Lifetime Cases 
show that for a wide range of lifetimes of H 2 0 2 in the soil, 
the atmospheric number density is affected only slightly. 
However, it is interesting to explore the behavior of the 
model for regimes where the soil may play a more active 
role in the destruction of H 2 0 2 . By identifying these re- 
gimes, it is possible to assess the effect of assumptions 
about the lower boundary condition used in more complex 
atmospheric models. In the loss regime case, the physical 
quantities are calculated for the condition that the soil 


and atmosphere are H 2 0 2 sinks of equal magnitude. This 
condition is met when the lifetime of H 2 0 2 in the soil is 

r = 8.8 hr. 

The resulting average surface concentration is very small, 
/7i = 0.275 nmoles cm 3 
and the extinction depth is very shallow 
z e = 3.5 x 10 3 cm. 

Thus we see that for the soil to represent a sink for H>0 2 
that is comparable with the atmosphere, surface concen- 
trations will be extremely small, with all the H 2 0 2 residing 
very close to the surface in a soil that efficiently catalyses 
its decomposition. 

4. DISCUSSION 

The nonlinear relationship between adsorbed and vapor 
phase H 2 0 2 results in a parabolic depth profile for [he 
solution to the soil diffusion equation. This curious result 
implies that H 2 0 2 is actually extinguished at a finite depth, 
rather than asymptotically approaching zero, as for a in- 
ear diffusion equation. Depth profiles for the Density and 
Lifetimes Cases are shown in Figs. 1 and 2, respectively. 
Figure 1 shows that when the model is constrained by an 
average surface concentration of 30 nmoles cm 3 , he 
h 2 o 2 is contained in a very shallow layer near the surfa;e, 
and no H 2 0 2 can penetrate below 1.5 cm. The lifeti Tie 
of H 2 0 2 in the soil is 4.3 terrestrial years, and the soil 
represents only about 1% of the total sink for H 2 0 2 . Very 
shallow extinction depths are, however, argued agaiist 
by the discovery of reactive soils under rocks (see Table 
I). It is likely, therefore, that the average surface concen- 
tration of H 2 0 2 is somewhat greater than the lower bound 
inferred from the Viking results. Figure 2 shows that when 
the lifetime of H 2 0 2 in the soil is constrained to be I0 5 
years, the average concentration of H 2 0 2 in the top 4 cm 
of soil is about eight times higher, at 237 nmoles cm' 3 . 
Again, the H 2 0 2 is confined to a shallow layer of 2.3 m 
near the surface. With this long a lifetime for H 2 0 2 in the 
soil, the soil becomes a negligible sink for H 2 0 2 . 

The soil diffusion equation is solved with a flux bound- 
ary condition at the soil-atmosphere interface, coupling 
it with the atmosphere mass-balance equation. In this 
way, given atmospheric production and loss rates, and soil 
loss and diffusion rates, concentrations may be predicted 
throughout the coupled soil-atmosphere system. Since 
the lifetime of H 2 0 2 against heterogenous catalytic de- 
struction under martian conditions is not well constrained , 
the model is characterized in terms of this parameter. 
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Lifetime (seconds) 

FIG. 3. The average concentration of H 2 0 2 in the top 4 cm of soil, 
plotted against the lifetime of H 2 0 2 . It can be seen that for lifetimes 
greater than 4.3 terrestrial years, the average concentration exceeds 30 
nmoles/cm\ 


Quantities of physical interest are the average concentra- 
tion of H 2 0 2 in the top 4 cm of soil, the extinction depth, 
and the atmospheric number density at the surface. These 
quantities are plotted as a function of the lifetime parame- 
ter in Figs. 3, 4, and 5. Figure 3 plots the average concen- 
tration of H 2 0 2 calculated to be in the top 4 cm as a 
function of the free parameter, the lifetime of H 2 0 2 in the 
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FIG. 4. The H 2 0 2 extinction depth plotted as a function ofthe lifetime 
of HiO z in the soil. The extinction of H 2 0 2 occurs at a finite depth due 
to the combined effects of adsorption and catalytic destruction retarding 
diffusion. For all but extremely long lifetimes (greater than 10 5 years), 
H 2 0 2 is unable to penetrate below 3 m. 



Lifetime (seconds) 

FIG. 5. Atmospheric concentration of H 2 0 2 plotted as function of 
the lifetime of H 2 0 2 in the soil. For lifetimes greater than 1 terrestrial 
day, soil losses do not substantially affect the concentration of atmo- 
spheric H 2 0 2 . As lifetimes in the soil become shorter, soil losses domi- 
nate the destruction of H 2 0 2 , and the overall atmospheric concentration 
is reduced. 


soil. Concentrations of 30 nmoles cm -3 , consistent with 
our interpretation of the Viking LR results, are achieved 
with lifetimes of 4.3 terrestrial years or more. The extinc- 
tion depth of H 2 0 2 as a function of the lifetime parameter 
is plotted in Fig. 4 . The relationship is nearly linear on 
a log-log plot, with depths to the unoxidized layer being 
less than 3 m for any lifetime less than about 10 5 years. 
In order to demonstrate the relative strengths of the atmo- 
sphere and soil as sinks for H 2 0 2 , the atmospheric number 
density is plotted as a function of the lifetime parameter 
in Fig. 5. The soil becomes a significant sink with respect 
to the atmosphere only if lifetimes are on the order of 
several hours or less. 

A significant result of the H 2 0 2 diffusion model pre- 
sented here is that for a very wide range of H 2 0 2 lifetimes 
in the soil, the oxidizing layer is confined to the top 3 m. 
This is due to the fact that as the partial pressure of H 2 0 2 
is reduced, the ratio of gas to adsorbate is also reduced. 
Therefore, as the mass flux drops, the role of adsorption 
in retarding diffusion becomes more significant. This re- 
sult is dependent upon the exact form of the adsorption 
isotherms that are incorporated into the model. Since 
laboratory data on the adsorption of H 2 0 2 under martian 
conditions is unavailable, the adaptation of laboratory- 
derived isotherms for water is clearly a weakness of the 
present model. However, simplified theoretical schemes 
of physical adsorption, namely Langmuir and BET iso- 
therms, roughly exhibit the same behavior — if vapor dif- 
fusivity falls off faster than partial pressure, the resulting 
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negative feedback will always limit H 2 0 2 diffusion to a 
fairly shallow depth. 

In calculating the results from the H 2 0 2 diffusion model, 
we have made several assumptions with regard to the 
behavior of H 2 0 2 in the martian soil. For the present 
analysis, the similarity in the dipole moments (12% differ- 
ence) and heats of fusion per gram (50% difference), for 
H 2 0 and H 2 0 2 serve to justify the adaptation of the labora- 
tory-derived isotherms. Although we believe that the ex- 
act nature of H 2 0 2 adsorption onto martian soil grains 
will not affect the basic results, a significantly better un- 
derstanding of the role that H 2 0 2 plays in the martian soil 
could be achieved through laboratory determinations of 
the appropriate adsorption isotherms. Similarly, labora- 
tory measurements of the stability of H 2 0 2 in Mars-like 
environments would be of great help in scaling the model 
to real Mars-like conditions. 

One possible explanation for the surprising stability of 
H 2 0 2 in the martian soil could be the observations that 
catalytic decomposition of peroxide is retarded at low 
temperatures. Solid H 2 0 2 is surprisingly inert. For exam- 
ple (Schumb et al. 1955), if added carefully so as not 
to allow thawing, chilled 0.5 N permanganate solution, 
particles of rust, and other catalysts were found to cause 
no decomposition of 90 wt% H 2 0 2 at -55°C. Decomposi- 
tion began only after melting started. The applicability of 
these results to adsorbed layers of H 2 0 2 is not clear but 
does suggest a possible reason for the low reaction rates. 

With a model of H 2 0 2 distribution in the martian soil 
and atmoshpere we may hope to address the following 
three questions: 

1 . Is a theoretical model of H 2 0 2 transport in the martian 
environment consistent with the Viking Biology experi- 
ment results? 

2. To what degree is the soil a sink for atmospheric 

h 2 o 2 ? 

3. At what depth may we expect to find an unoxidized 
layer, possibly containing organic compounds? 

In order for the model to be consistent with some inter- 
pretation of the Viking Labeled Release experiment, we 
may constrain the average concentration of FI 2 0 2 in the 
top 4 cm of the soil to be at least 30 nmoles cm" 3 . This 
may be achieved in the model if the lifetime of H 2 0 2 in 
the soil is at least 4.3 terrestrial years. In this sense, the 
model is not inconsistent with the additional inference of 
a H 2 0 2 lifetime of 10* years. Additionally, at these soil 
concentrations and lifetimes, we see that the soil is a poor 
sink for H 2 0 2 . Therefore, atmospheric concentrations of 
H 2 0 2 predicted by this model are consistent with more 
complex atmospheric models that assume a passive sur- 
face. We may conclude that as long as the lifetime of 
H 2 0 2 in the soil is greater than 4.3 terrestrial years, the 
model is consistent with Viking Biology experiment re- 


sults. For any value of lifetime of H 2 0 2 in the soil tha is 
consistent with the Viking Biology experiment results, he 
soil is a poor sink for H 2 0 2 compared with the atmosphere. 
Unless the lifetime is on the order of several hours or 
less, the atmosphere may be considered the primary s nk 
for H 2 0 2 . Given this condition, more complex models of 
H 2 0 2 production and loss in the martian atmosphere fiat 
neglect losses to the soil do so without loss of accuracy. 
For all lifetimes shorter than 10* years, the depth to com- 
plete extinction of H 2 0 2 was less than 2.3 m. Although it 
might be expected that higher diffusion depths can be 
achieved by a gas, adsorption and subsequent catalytic 
destruction severely limit the depth to which H 2 0 2 may 
penetrate. As Chyba et al. (1989) have pointed out, small- 
scale cratering has probably impact gardened the regolith 
to a depth of 5 m or more. Therefore, if H 2 0 2 is he 
primary agent responsible for the oxidation of organics 
in the martian regolith, it may be inferred that the depth 
at which organics may still exist is not limited by he 
existence of H 2 0 2 . It is possible that organic material may 
remain intact beneath a shallow impact-gardened layer. 

Recently, Wright et al. (1989) have analyzed the mete- 
orite EETA 79001, one of the SNC meteorites believed 
to have originated on Mars (Becker and Pepin 1984, Bo- 
gard et al. 1984), and found evidence that it contains more 
than 400 ppm degradable organic material. The isotopic 
ratio of carbon in the organic material is suggestive of 
a biological origin, although terrestrial contamination in 
these measurements cannot be ruled out (Jakosky 19S t). 
The organic material may be derived from organic bearing 
sediments that were laid down during an earlier, more 
clement period when life was widespread on the martian 
surface (McKay and Stoker 1989). Since SNC meteorites 
are thought to have resided on Mars up to 180 mill on 
years ago, and conditions on Mars that were conducive 
to life had disappeared by 2 billion years ago (Pollack et 
al. 1987, McKay and Stoker 1989), it is possible that he 
organic material seen in the meteorite was formed during 
an earlier epoch. Using the proposed mechanism (Melosh 
1985) to explain how the SNC parent body was ejected 
from Mars, one can estimate the approximate depth ie- 
neath the martian surface from which the meteorite ETA 
79001 originated. Because the meteorite contains organ- 
ics, it presumably originated from below the oxidizing 
soil layer. The original depth of the material is estima ed 
from the fact that the meteorite was not exposed to cosmic 
rays in space and therefore must have been shielded by 
a minimum of 3 m of material (Bogard etal. 1984). Because 
the meteorite material is only weakly shocked, it mast 
have come from near the martian surface. The maxim jm 
depth of material ejected from Mars depends on the size 
of the impacting object that ejected it (Melosh 1985). For 
an impactor of 10 km, the maximum depth of the material 
ejected is approximately 40 m (Melosh 1985). Thus pre- 
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suming that the SNC meteorites are from Mars, unoxi- 
dized organic carbon should be found at reasonably shal- 
low depths, a few tens of meters at most. This is in 
agreement with the results calculated from our soil diffu- 
sion model. 

For the purpose of discussing the oxidizing soil on 
Mars, we have adopted a “standard” model of martian 
soil oxidants. The results of the Viking biology experi- 
ments are assumed to be explained by the existence of 
three oxidants — a strong, thermally stable oxidant or ad- 
sorbed 0 2 is invoked to explain the GEx results, H 2 0 2 is 
assumed to be responsible for the rapid decomposition of 
organics in the LR experiment, and an additional weak 
oxidant is required to explain the slow continued release 
of C0 2 in the GEx and LR experiments (Klein 1978, 1979). 
The purpose of the diffusion model presented here is to 
explore the transport of one of these oxidants, H 2 0 2 and to 
see to what extent it may affect the depth to an unoxidized 
layer. Given the assumptions of the model and the rough 
extrapolations from laboratory data, we find that for H 2 0 2 
lifetimes of less than 10 5 years, the depth to the unoxidized 
layer is 3 m or less. 
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